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On mode conversion and wave reflection in magnetic Ap 
stars 



00 

o 
o 

(N 
D 



O 



> 

m 

00 

i> 
O 

(N 
O 
oo 
O 



S. G. and M. S. Cunha^ 

^ Centra de Astrofisica da Universidade do Porto, Rua das Estrelas, s/n, 4150-762 Porto, Portugal 
^ Departamento de Matemdtica Aplicada, Faculdade de Ciencias da Universidade do Porto, Portugal 



Accepted ??. Received ??; in original form 2007 July 



ABSTRACT 

We investigate the effect of a strong large scale magnetic field on the reflection of high 
frequency acoustic modes in rapidly oscillating Ap stars. To that end we consider a toy 
model composed of an isothermal atmosphere matched onto a polytropic interior and 
determine the numerical solution to the set of ideal magneto-hydrodynamic equations 
in a local plane-parallel approximation with constant gravity. Using the numerical 
solution in combination with approximate analytical solutions that are valid in the 
limits where the magnetic and acoustic components are decoupled, we calculate the 
relative fraction of energy flux that is carried away in each oscillation cycle by run- 
ning acoustic waves in the atmosphere and running magnetic waves in the interior. 
For oscillation frequencies above the acoustic cutoff we show that most energy losses 
associated with the presence of running waves occur in regions where the magnetic 
field is close to vertical. Moreover, by considering the depth dependence of the energy 
associated with the magnetic component of the wave in the atmosphere we show that 
a fraction of the wave energy is kept in the oscillation every cycle. For frequencies 
above the acoustic cutoff frequency such energy is concentrated in regions where the 
magnetic field is significantly inclined in relation to the local vertical. Even though our 
calculations were aimed at studying oscillations with frequencies above the acoustic 
cutoff frequency, based on our results we discuss what results may be expected for 
oscillations of lower frequency. 

Key words: Stars: oscillations - Stars: variables - Stars: magnetic. 



1 INTRODUCTION 

A number of rapidly os cillatinR Ap stars , including the 
two prototypes a Cir l|Kurtz et al.l and HR1217 

l|Kurtz et al.l l2005l ). are known to pulsate with frequencies 
that exceed the acoustic cutoff frequency expected from ap- 
propriate stellar models. This fa ct has been a matter of 
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long debate over the years (e.g. 



Shibahashi fc Saio 



1998h . 



lAudard et all Il998l : ICunhal 11998 : iGautschv et al.l 
Generally, it has been suggested that the apparent dilemma 
posed by the observation of oscillations of such high fre- 
quency in roAp stars results from the use of a Tempera- 
ture - Optical depth (T — r) relation that is inadequate for 
these peculiar stars. However, except for the case in which 
the T — T relation was modified to simulate the presence 
of a chromosphere, for which no observational evidence ex- 
ists, models with modified T — t relations considered in the 
above mentioned works, have not been able to fully resolve 
the apparent discrepancy. 
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Even though the oscillations observed in roAp stars 
have essentially an acoustic nature, it is nowadays well 
established that the latter acquire a magnetoacoustic na- 
ture in the outer layers of the star, where the magnetic 
and gas pressure become comparable. A number of stud- 
ies investigating on the effect of the magnetic field on 
the properties of the oscillations, such as eigenfrequen- 
cies and ei genf unctions, were carried out by different au- 



thors (e.g. IPzie mbows ki fc Good3 199(: ^ Bigot et al 



iCunha 



.20001 : 
20061 ). 



ICunha fc Goughll2000l : ISaio fc Gautschvir2004 

Mode conversion - i.e, the process by which the acoustic 
wave in the interior passes part of its energy to waves of a 
different nature through the coupling with the magnetic field 
in the region where the gas and magnetic pressures are com- 
parable - was considered in all the above mentioned works. 
In some of them the waves were artificially reflected at the 
surface of the star. In those cases only the energy passed 
onto running magnetic waves in the interior was considered. 
In the studies of Dzicmbows ki & Goodc (1996) and Cunhal 
(|2006l ). on the other hand, the authors did not make the 
former assumption, and as a result, part of the energy asso- 



2 S. G. Sousa and M. S. Cunha 



ciated with the acoustic wave in the interior was allowed to 
be converted, through the coupling with the magnetic field, 
into running acoustic waves in the high atmosphere, as well 
as into running magnetic waves in the interior. 

While incorporating the process of mode conversion, the 
works mentioned above aimed mostly at studying the overall 
effect of the presence of the magnetic field on the eigenfre- 
quencies and eigenfunctions of the oscillations. Thus, in no 
case the impact that such process might have on the issue 
of mode refiection at the surface of roAp stars was prop- 
erly analysed. This issue, which is directly related to the 
long debated problem of the cutoff frequency in roAp stars 
discussed at the beginning of this Section, is studied in the 
present work in a toy model composed of an isothermal at- 
mosphere matched onto an index 3 polytropic interior. We 
note that the phenomena of mode conversion is known to 
take place in other classes of pulsating stars. In particular, 
in recent years a number of studies related to mode conver- 
sion were carried out also in the contex t of solar pulsations 
in reg ions of strong magnetic field (e.g. iMcDougall fc HoodI 
I2OO7I . and references therein). 
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Figure 1. Schematic view of the the outer layers of the magnetic 
star. The two outer regions, where the magnetic pressure is larger 
than, or of the same order of the gas pressure will be denominated 
magnetic boundary layer. 



2 THE MAGNETIC BOUNDARY LAYER 

2.1 Governing equations and general assumptions 

In this work we are particularly concerned with pulsations in 
Ap stars. Observations indicate that Ap stars have intense, 
large scale, ma gnetic fields, often with predominant dipola r 
structure fe.g. iHubrig et al.ll2005l . 12004 : IWade et al.ll2000h . 
Accordingly, in the present work we will consider that our 
stellar model is permeated by a dipolar magnetic field. Since 
this is a force free field, it does not influence the equilibrium 
structure of the star, which therefore is governed by the 
system of equations. 



— - = • Vpo - PoV • v, 

ot 



dt ' 
dt ' 



(1) 



(2) 



(3) 



where po is the density, po is the pressure, go is the grav- 
itational field, and the subscript refers to the equilib- 
rium quantities. In the above equations d stands for partial 
derivatives and t is the time. 

In the limit of perfect conductivity, adiabatic motions 
are governed by the set o f ideal mag neto-hydrodynamic 
(MHD) equations (see, e.g.. |Priestlll982l ). Considering small 
deviations from the equilibrium state and then linearizing 
the MHD equations by neglecting squares and products of 
the small quantities (that will be denoted with the subscript 
1), and neglecting the Eulerian perturbation to the gravita- 
tional potential (Cowling approximation), we arrive to the 
following system of equations. 



dv 
'dt 



^ = Vpi + ^Vpo + (V X Bi) X Bo, 



po 



Po 



ppo 



dpi „ „ 

— — = -V ■ Vpo - 7PoV ■ V, 
ot 



(4) 



(5) 



dt 



= V X (« X Bo) 



V-Bi =0, 



(6) 

(7) 
(8) 



where v is the pulsation velocity, 7 is the the first adia- 
batic exponent, is the magnetic permeability, and B is 
the magnetic field. The perturbed Lorentz force in equation 
Q expresses the magnetic field influence on the oscillations. 

In our study we will deal only with the outer layers of 
the star. We shall consider three different regions (cf. Fig.[T]), 
namely, 

• A region in the atmosphere where the magnetic pressure 
dominates over the gas pressure {/3~ region); 

• A region where the magnetic and the gas pressure are 
of the same order of magnitude; 

• A region where the magnetic pressure is negligible when 
compared with the gas pressure {f}^ region). 

Here /3 is defined as the ratio between the gas pressure, pg, 
and the magnetic pressure, pm ~ 8^/2/1. The region f3~ 
corresponds to the layers where /3 <C 1, while the region /3"^ 
corresponds to the layers where /3 ^ 1. 

We shall refer to the two outer regions as the magnetic 
boundary layer. We note that since in the interior of the star 
the gradient of the gas pressure is much greater than the 
Lorentz stresses, the dynamics there is effectively field free. 
Thus, the magnetic field influences the global oscillations 
only through its direct contribution to the restoring force in 
the magnetic boundary layer. 

In this study we will consider perturbations to the equi- 
librium state of the model associated with high order, low 
degree, axisymmetric pulsations only. These are the pulsa- 
tions that are typically observed in roAp stars. Moreover, 
we will not consider pure incompressible perturbations (even 
though in the P'^ region the decoupled magnetic component 
becomes almost incompressible). 
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Since the magnetic field is assumed to have a dipolar 
configuration and the pulsations are axisymmetric, the per- 
turbations depend on co-latitude 9 but not on the axial angle 
0. 

Under these conditions the system of equations (j4])-(|8]) 
admits solutions of the form, 



(9) 



where ^ is the displacement vector defined hy v — lo is 
the oscillation frequency and r is the radial coordinate, in 
the spherical coordinate system (r, 9, <j)). 

Since the boundary layer is thi n and the magnetic field 
varies only on large scales, we follow iDziembowski fc Goodd 
l|l996ll and ICunha fc Goughl (|200(]| ) and solve our problem 
locally by adopting, at each latitude, a plane parallel ap- 
proximation, and assuming a locally uniform field. 

At each given latitude, our Cartesian system is defined 
such that the local vertical coordinate z increases outwardly 
and is zero at the surface of the star. The coordinates in the 
horizontal plane, x and y, are chosen such that the x axis is 
parallel to the horizontal component of the magnetic field. 
Hence Bo = (B.,0,B,). 

In accordance with our approximation of a locally uni- 
form magnetic field, we disregard the derivatives of Bx and 
-Bz in equations (|4ll-([8]). Since the oscillations considered 
are axisymmetric, d^/dy = 0. Moreover, given that the ob- 
served modes are of low degree, we expect that ^ will vary 
much more rapidly in the z direction than in the x direc- 
tion. Hence, we also approximate d^/dx « 0. The latter 
approximation is exact for pulsations that in the absence of 
the magnetic field would be spherically symmetric (i.e. with 
degree 1 = 0). 

The presence of a magnetic field, and the consequent 
Lorentz force that is generated, introduces horizontal motion 
in the magnetic boundary layer even in the case of modes 
that in the absence of a magnetic field would be radial. We 
can decompose the displacement into its vertical component 
and its horizontal component defined, respectively, by = 
^ • fiz and (^x = ^ • e^,, where and are unit vectors 
in the z and x directions, respectively. The fact that the 
horizontal displacement is along is a consequence of the 
axisymmetry of the problem. 

Using equations ©-((Tl) to eliminate pi, pi and Bi from 
equation (Q, neglecting the derivatives of Bx and B^, and 
finally separating into vertical and horizontal components 
we obtain. 



■ uj''po^. = [(Bo ■ V)(V ■ i)] + -{Bo ■ Vfix, (10) 



-^W. = ^- — [(So-V)(V-0] 



dz /i 
fi po dz 



(11) 



where. 



W = i-Vpo + (7P0 ^)V • i - -{Bo ■ V)(Bo • ^). (12) 
M P 



2.2 Dimensionless Variables 

When solving our problem we will make the system of equa- 
tions pOII - Hlip dimensionless by defining the new variables: 



ri = -z/R, 
p = po/po, 



- Po/po, 



£x — £,x/R, Ez — (.z/R, 



bi = Bi/{fipouioR' 



2 R2^1/2 



C2 = po/{po<^oR^) 



where po is the pressure at the stellar surface, po the den- 
sity at the stellar surface and ujo = \/GM/R? , where M 
and R are, respectively, the mass and radius of the star, 
and G is the gravitational constant. Following the same 
notation, the dimensionless magnetic field vector becomes 
''o = (fei,0,6z) and the total dimensionless displacement 
becomes e = Sx&x + EzEz. Moreover, we define the dimen- 
sionless sound speed as c = {C2'yp/ pY^"^ ■ 

Combining these new variables with the system of equa- 
tions pop - (|lip and simplifying, we obtain: 



hxhzS^ + h^Zx + o" P^x 



0, 



(13) 



(7C2P + 6^)e'z' + 'yC2pe'^ + a'^pe^ - bxh^e" = 0, (14) 

where the prime ' stands for derivatives in order to the new 
independent variable rj. 



2.3 Magnetic coordinate system 

When the system of equations (|13p - (|14p is studied in the 
extreme case of no gas pressure it is found that the displace- 
ment is perpendicular to the magnetic field. This is a direct 
consequence of the Lorentz force being the only restoring 
force playing a role in that extreme case. With this in mind 
we will define a second coordinate system with axes directed 
along and perpendicular to the magnetic field (see Fig. 
To change from the first to the second coordinate system we 
use the relation: 



£x = e|| cos a — ei_ sin a 
Ez = £_L cos Q -I- Ell sin a 



(15) 



where a is the inclination angle of the magnetic field in re- 
lation to the horizontal direction and and ex are the di- 
mensionless components of the displacement along and per- 
pendicular to the magnetic field direction, respectively. This 
coordinate system will be useful when dealing with the prob- 
lem in the /3~ region. On the other hand it will be useful to 
use the {x, z) coordinate system in the region, since in 
the region where the gas pressure dominates over the mag- 
netic pressure the vertical compo nent of the di splacement 
is essentially an acoustic wav e (e.g. ICampbell fc Papaloizoul 
1 19861 : ICunha fc Goughll2000l ) . 



3 STELLAR MODEL 

The main goal of this work is to understand how high fre- 
quency magnetoacoustic pulsations are reflected in the outer 
layers of Ap stars, and what fraction of the pulsation energy 
is lost through running waves. Thus, we start by consid- 
ering a simple model that still contains the main physical 
ingredients necessary to study this problem, namely a plane 
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Figure 2. Schematic view of the components of the displacement 
e projected onto the two coordinate systems used in the present 
work, ex and Sz are the dimensionless horizontal and vertical 
components of the displacement, respectively, ey and e±_ are the 
dimensionless components of the displacement along and perpen- 
dicular to 60 ) respectively. The angle between the local magnetic 
field direction and the local horizontal coordinate x is designated 
by a. 



parallel envelope under constant gravitational acceleration, 
composed of an isothermal atmosphere matched at the sur- 
face onto a polytropic interior of index 3. We hope to refine 
our analysis using a more sophisticated stellar model, in fu- 
ture work. 

All results presented are for a model of mass M = 
2.0 Mq and radius R = 2.0 R©, where R0 and M© are 
the radius and mass of the sun, respectively. The temper- 
ature of the isothermal atmosphere (where the perfect gas 
law is assumed to hold) is T = 8440 K and the atmosphere 
is matched continuously onto the polytropic interior at a 
pressure P = 7563 gcm~^ . These values were chosen by 
comparison with a main sequence (age = 468Myrs) stellar 
model of th e same mass and radius computed with the CE- 
SAM code (Morel" igQ?"). 

The expr essions for th e density and pressure profiles are 
given by (e.g. lGoughlll99 
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Figure 3. Characteristic dimensionless magnetic (dashed line) 
and acoustic (full line) wavelengths as a function of dimensionless 
depth in the star, as defined by equations 12011 and (I2H . respec- 
tively. These were derived for the model under study, considering 
a cyclic frequency v = 3.08 mHz and a characteristic magnetic 
field of 3000 G. 



under study defined, respectively, by 

27ri6ol 



and 



A. 



(7^ 
2-KC 

a 



(20) 



(21) 



assuming a cyclic oscillation frequency v = 3.08 mHz and 
a magnetic field strength of 3000 G. In the isothermal at- 
mosphere [j] < 0) the characteristic magnetic wavelength is 
much greater than the characteristic acoustic wavelength. 
Deeper in the star there is a region where the two wave- 
lengths are of the same order of magnitude, and deeper still 
a region where the acoustic wavelength is much greater than 
the magnetic wavelength. Thus, it is expected that the mag- 
netoacoustic mode described by equations (|10|l - pi|l will de- 
couple onto acoustic and magnetic components, both suffi- 
ciently high in the atmosphere and sufficiently deep in the 
interior of the star, so far as the waves propagate there. The 
exact regions of decoupling will depend on the magnetic field 
intensity and on the gas pressure in the outer layers of the 
star. 



:e«. 



(16) 



p = e«, (17) 

in the isothermal atmosphere, where Ti. = H/R and H is the 
scale height given hy H — poR^ /{poGM), and by 

P^(j + l)^ (18) 

P^(j + l)\ (19) 

in the polytropic interior, where rjs = 4Ti.. 

In Fig. [3] we show the characteristic (dimensionless) 
magnetic and acoustic oscillation wavelengths for the model 



4 DECOUPLING REGIONS 
4.1 /3+ Region 

As shown by [Roberts fc SowardI l|l983D . deep enough in the 
star the magnetoacoustic wave decouples into its fast acous- 
tic and slow magnetic components, and it is possible to find 
an approximate solution for the horizontal component of 
the dimensionless displacement Ex, by applying appropriate 
standard asymptotic theory valid for high frequency modes. 
Within the approximations considered in the present work, 
the horizontal component of the displacement in the re- 
gion is entirely magnetic and has the form. 
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and 



(22) 



where D and E are depth independent amphtudes. Follow- 
ing the same authors we assume that this component of 
the displacement is an inwardly propagating wave and set 
E = 0. This is justified by the fact that the latter describes 
a magnetic wave that is expected to dissipate while propa- 
gating inwardly, due to the rapid increase of its wavenumber 
with depth. 



4.2 f3~ Region 

If the oscillation frequency is sufficiently high, the wave will 
propagate in the atmosphere. In that case we expect the dis- 
placement to decouple into a slow acoustic and a fast mag- 
netic component, as a result of the difference in the char- 
acteristic magnetic and acoustic wavelengths in this region. 
Thus, under these conditions the total displacement may be 
written as the sum of fast (subscript /) and slow (subscript 
s) oscillatory components. 



(23) 



In the l3~ region the perturbed Lorentz force dominates the 
restoring force for perturbations that are perpendicular to 
the unperturbed magnetic field. Hence, the component of 
the displacement in the direction perpendicular to the un- 
perturbed magnetic field will be associated essentially with 
a compressible magnetic wave. However, to first order the 
Lorentz force has no effect for perturbations along the di- 
rection of the unperturbed magnetic field. Thus, the dimen- 
sionless displacement along that direction will be associated 
with a wave that is essentially acoustic. We thus expect to 
have, 



and 



> ex. 



(24) 



(25) 



Rewriting the system of equations (|13p - p4|l in the co- 
ordinate system defined by equations (|15p we find. 



7C2p(e'I+eI|'^^ )+■ 



I sm a \ 
+ei +- 



0, 



sm a cos a 



e" + (7^p^£|| cos a - 



0. 



(26) 
(27) 



Since the coefficients multiplying the displacement in equa- 
tions (|26p - (|27p are not wave like functions, these equations 
must be satisfied also by the fast and slow components of 
the displacement separately. Neglecting exs when compared 
with e\\g and e\\f when compared with ex/, we find that the 
equations governing the slow and fast components of the 
displacement are, within the referred approximations. 



+ ■ 



+ 



^C2P sin a 



■ 



(28) 



respectively, where 



+ ■ 



7C2 cos^ a[p+ 13'^) 



ex/ = 0, (29) 



7C2 cos^ a 

This decoupling of the magnetoacoustic wave into acoustic 
and magnetic components in the P~ region will be verified 
numerically in Section [S] Equations (|28p and (|29|) may be 
written in the form of a standard wave equation. To that 
end we define the new variable 



1/2 
= P ^11 



and introduce it in equation (|28|) to find, 
where 



P_\ 

2p J 2 p 7C2P sin^ a 



(30) 
(31) 

(32) 



For equation (|29|) we can use a similar variable, defined by 
Hx/ = (P + B')''"ex/, (33) 



to find, 



where 



-X/ + «x- 



x=x/ 



p 



0, 



(34) 



.2(p + B2)/ 2(p + B2) jC2{p + B'^) cos'^ a 

(35) 

Equations (|28p and (|29|) were derived under the assump- 
tion that two decoupled wave components with significantly 
different characteristic wavelengths can be identified in the 
/3~ region. That condition is satisfied when fcy and k'^ are 
both positives, fcy is positive for oscillations with frequencies 
larger than CTc where, 

'^c = ^sinQ. (36) 

Except for the dependence in the angle a, this expression 
is similar to the cutoff frequency for acoustic waves in an 
isothermal atmosphere with no magnetic field (hereafter 
named acoustic cutoff frequency) . In fact, when the magnetic 
field is vertical the local vertical motion associated with the 
radial pulsations (or a low degree mode) will not be affected 
and we recover the acoustic cutoff frequency. As the mag- 
netic field becomes more inclined in relation to the local 
vertical, the necessary frequency for propagation of acoustic 
waves in the atmosphere decreases. 

Concerning the magnetic component we find that k\ 
is positive for oscillations with frequencies larger than am, 
where 

-| 1/2 

1 1 

2 



n 



4 1 



(37) 



The frequency am increases with height in the atmosphere, 
taking a maximum value (as ^ — oo) of 



V2c 

o-m = cos a. 



(38) 
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Figure 4. The acoustic (full line) and magnetic (dashed line) fre- 
quencies defined by equations II36I I and II38I I. respectively, for our 
model, and their dependence on the inclination of the magnetic 
field in relation to the local horizontal direction, a. 



In Fig. 13] we show the frequencies ctc and (T„ifor the model 
used in the present work. 

We emphasize that equations (|28p and (|29p are approxi- 
mately valid only in the outer layers, where the magnetic and 
acoustic components of the wave are decoupled. Hence, they 
may not be used to directly infer under which conditions the 
magnetoacoustic wave is reflected (or partially reflected) in 
the layers below. For that reason we avoided using the term 
'critical frequencies' when referring to ctc and Om- In sec- 
tion |S] we shall look at the problem of wave reflection and 
pulsation energy losses more closely, using the numerical so- 
lutions obtained for the magnetic boundary layer. 

The components of the displacement in the /3~ region 
can be obtained from the solutions to equations (|28|) and 
\29\ . The general solution to equation (|3H) is given by 



(39) 



where A and B are depth independent amplitudes. Since no 
energy can be sent in from outside the star, this must repre- 
sent a wave propagating outwardly through the isothermal 
atmosphere. Thus B must be zero. Using equation (|30|) we 
thus find the solution to equation (|28p to be. 



P 



A 



ikuT) 



(40) 



To find an approximate solution to equation (I29|l we 
first note that when terms of order are neglected, fc^ 
becomes 



fc 



|boP 



- 2 
0"m 



(41) 



Using p as the independent variable, equation (I34|l thus be- 
comes. 



ad^Si^ dHx/ 



dp2 



+ P- 



dp 



7f 



(a^ - &l) E^f = 0. (42) 



The solution to equation (|42p c an be written in terms 
of Be ssel functions Jo and Yo (e.g. lAbramowitz fc StegunI 
Il972l . pg 362), namely. 



CJo (2^pj + DY„ , 



(43) 



where C and D are depth independent amplitudes and X = 
7i^|6o| (f^ — (Tm). Since Yo diverges as p 0, -D must 
be zero. Thus, using equation (|33p . we find that within the 
approximations described, the solution to equation (|29p has 
the form, 



C 



(44) 



5 SOLUTIONS IN THE MAGNETIC 
BOUNDARY LAYER 

Throughout most of the magnetic boundary layer, the dis- 
placement is associated with a magnetoacoustic wave which 
cannot be described, not even approximately, by decoupled 
magnetic and acoustic components. Thus, to determine the 
form of the displacement in this case it is necessary to solve 
the system of equations p3|l - (|14| l numerically applying, si- 
multaneously, appropriate boundary conditions. In the ab- 
sence of a magnetic field, the acoustic cutoff frequency is 
equal to the polar value of CTc- Since our goal is to investi- 
gate on the reflection of oscillations of which frequency is 
above this acoustic cutoff frequency, in what follows we will 
consider the case of a cyclic frequency well above c/(2n2Ti.). 



5.1 Boundary Conditions 

Since the system of equations to be solved is linear, in order 
to solve it we need to apply three boundary conditions. Two 
of these conditions will be applied in the (3~ region while 
the third will be applied in the region. 

At the outermost point of the atmosphere we match 
the perturbed magnetic field into a vacuum fiel d. This con- 
dition implies that e'± — (e.g. ICunhal [20061 ). To derive 
the two additional boundary conditions we use the fact that 
two approximate analytical solutions can be obtained, one in 
each of the decoupling regions. In the /3~ region the parallel 
displacement, ey, is essentially associated with the acoustic 
wave. Hence, in this region the parallel component of the 
numerical solution must be well described by the analyti- 
cal expression given by equation (|40|l . Moreover, in the /3+ 
region the horizontal component of solution is essentially as- 
sociated with the magnetic wave. Hence in this region the 
horizontal component of the numerical solution must be well 
described by the asymptotic solution given by equation (|22|l 
with E ^0. 

To find a solution that obeys both the system of equa- 
tions p3p -(|14 p and the boundary conditions given above we 
start by finding three linearly independent particular solu- 
tions that satisfy the system of equations and the first of 
these boundary conditions. To generate these three partic- 
ular solutions we had to provide starting conditions that 
assured they were linearly independent. The first, second 
and third solutions were generated by assuming at the out- 
ermost point of integration, respectively, {e± = 1, ey — 0, 
ejl = 0), (e± = 0, £|| = 1, eJi = 0) and (e± = 0, £|| = 0, 
£|l = 1). To apply the two additional boundary conditions 
we write the total solution as a linear combination of these 
particular solutions. 
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£ = si + cie2 + C2es (45) 

The complex coefficients ci and C2 are determined by 
applying the two remaining boundary conditions, i.e. by 
matching simultaneously the relevant components of the nu- 
merical solution to the analytical solution (|22p with E — 
in the deeper layers and to the analytical solution (|40|) in 
the atmosphere. When performing these matchings we in- 
troduced an integer index g which reflects the positions at 
which the matchings are carried out in a grid of depths. 
The value of this index decreases from a thousand to zero 
as the depths at which the matchings are carried out are 
moved from the region where /3 ~ 1 to the regions where 
the magnetic and acoustic wave components are decoupled 
(see Appendix A for details). As mentioned in the previous 
Section, we can check the decoupling of the magnetoacoustic 
wave in the /3~ and /3"^ regions numerically. If such decou- 
pling occurs, the coefficients ci and C2 must tend to constant 
values as the index q tends to zero. 

Fig.[S]shows the behaviour of ci and C2 when the match- 
ings are carried out at different depths. The results shown 
are for an oscillation with cyclic frequency v — 3.08 mHz 
and a dipolar magnetic field of polar strength Bp — 3000 G 
at a latitude such that its inclination with respect to the 
local horizontal axis is a = 45°. When the matchings are 
carried out in the region where /3 ~ 1, which in Fig. [5] cor- 
responds approximately to q ranging between 600 and 1000, 
the coefficients ci and C2 vary significantly with the index q. 
This is because the analytical solutions used in the match- 
ings are not good representations of the true solutions in 
the regions considered. However, as we move the matchings 
simultaneously to the regions /3~ and (close to q — 0) 
ci and C2 tend to constant values, confirming the expected 
decoupling of the magnetoacoustic wave in both regions. 

5.2 Numerical Solutions 

Having calculated the constants for the linear combination 
we can determine the total solution for the displacement. For 
illustration we show in Fig. |6] the real part of dimensionless 
displacement for a dipolar magnetic field of polar strength 
Bp = 3000 G, at a latitude such that a = 45°, and a cyclic 
oscillation frequency v = 3.08 mHz. The figure shows the 
components of the dimensionless displacement perpendicu- 
lar and parallel to the direction of the magnetic field (upper 
panels) and the components of the dimensionless displace- 
ment in the vertical and horizontal directions (lower panels). 

The decoupling of the magnetoacoustic wave into acous- 
tic and magnetic components both in the (3^ and /9~ regions 
can be seen in the plots. To illustrate that, zooms of Fig.|S]in 
the interior and in the atmosphere are shown in Fig. [7] and 
Fig. [S] respectively. Fig. [7] shows the vertical (upper panel) 
and horizontal (lower panel) solutions in the interior. The 
decoupled acoustic and magnetic waves are readily seen. The 
horizontal component is a rapidly varying magnetic oscilla- 
tion, while the vertical component is an acoustic oscillation 
with a significantly larger wavelength. Moreover, inspection 
of the real and imaginary parts of the horizontal displace- 
ment confirm the running nature of the magnetic wave in the 
interior, characterized by real and imaginary parts of simi- 
lar amplitude and out of phase by 7r/2. The nearly standing 



nature of the acoustic wave in the interior is also clear from 
the imaginary part of the vertical displacement, which am- 
plitude tends to zero with increasing depth (alternatively, 
a standing wave could have real and imaginary parts with 
significant amplitude, but in phase). Fig. |8] shows the com- 
ponent of the dimensionless displacement perpendicular to 
the magnetic field direction (upper panel) and the compo- 
nent of the same vector along the magnetic field direction 
(lower panel), in the atmosphere. As before, the acoustic 
and magnetic components are easily identified. The compo- 
nent along the direction perpendicular to the magnetic field 
tends to a constant, as expected for the magnetic wave from 
the analysis carried out in Section [4. 2 1 while the component 
of the dimensionless displacement along the direction of the 
magnetic field is an acoustic oscillation with a wavelength 
comparable to the thickness of the atmospheric layer shown. 
Once again, inspection of the real and imaginary parts of 
the displacement confirm the running nature of the acous- 
tic wave in the atmosphere, and the standing nature of the 
magnetic wave in the same region. We note that although 
the standing nature of the magnetic wave in the atmosphere 
is connected to the condition of matching of the perturbed 
magnetic field onto a vacuum field, the form of this compo- 
nent of the displacement is rather insensitive to the place at 
which that boundary condition is applied, so far as the latter 
is applied sufficiently high in the atmosphere. Also, as shall 
be discussed in Section|Sl this nature of the magnetic wave is 
not associated to a full wave reflection at a particular place 
in the atmosphere, but rather to a progressive decrease of 
its energy content with height. 

Fig. [S] [3 and |S] illustrate also the fact that in the at- 
mosphere the decoupling takes place in the parallel and per- 
pendicular directions, while in the interior it takes place in 
the horizontal and vertical directions, instead. 



6 ENERGY FLUXES 

6.1 Flux carried by running waves 

The acoustic waves that propagate outwardly in the atmo- 
spheres of roAp stars and the magnetic waves that propa- 
gate inwardly in their interiors, both take energy away from 
the global oscillations. This is true both for oscillations with 
frequencies below and above the acoustic cutoff frequency. 
For oscillations with frequencies below the acoustic cutoff 
frequency, energy losses through running waves in the at- 
mosphere are restricted to regions where the inclination of 
the magnetic field is such that u > a^- On the other hand, 
oscillations with frequencies above the acoustic cutoff are 
expected to lose energy through running waves in the at- 
mosphere regardless of the inclination of the magnetic field. 
For this reason the observation in some roAp stars of oscil- 
lations with frequencies above the acoustic cutoff frequency 
has been a matter of long debate. At first sight, one would be 
led to think that perturbations of such high frequency would 
simply propagate away in the atmosphere and be dissipated 
before they could grow to observable amplitudes. However, 
the simple fact that a global oscillation of such high fre- 
quency is observed indicates that only part of the mode en- 
ergy is removed by the running acoustic and magnetic waves. 
Moreover, it also indicates that the energy input to such a 
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Figure 5. The figure shows the behaviour of the complex coefliicients of the linear combination expressed by equation II45II . A high value 
of the index q means that the matchings are carried out in a region where /3 ~ 1, while low values of that index mean that the matchings 
are carried out in the regions /3+ and (3~ . ci and C2 tend to constant values as the index q tends to zero, confirming the decoupling of 
the magnetoacoustic components both in the outer atmosphere and in the interior. Results are for a dipolar magnetic field with polar 
strength Bp = 3000 G at a latitude such that its inclination with respect to the local horizontal axis is a = 45°, and a cyclic frequency 
V = 3.08 mHz. 



mode in each cycle, through the opacity mechanism acting 
in th e region where hydrogen is ionized l|Balmforth et al] 
I2OOII : [Cunha 2002 : Saio .200 5). is sufficient to compensate 
all sources of mode energy losses, including the energy re- 
moved by the running waves. 

With the question of mode reflection in mind, in this 
Section we estimate the fraction of energy flux that is car- 
ried away by running waves and its dependence on the in- 
clination of the magnetic field. As before, the calculations 
are carried out in a local plane-parallel approach, assuming 
a dipolar magnetic field and using the numerical solutions 
for the displacement obtained in Section [5] to compute the 
relevant energy fluxes at each co-latitude. 

Let us consider a wave that propagates outwardly from 
the interior of the star carrying a given energy flux (see 
Fig. [9]). At some point the wave is reflected, or partially re- 
flected, and propagates inwardly again, back to the starting 
point. As a result of the magnetoacoustic coupling near the 
surface layers, this wave transfers part of its energy to a 
inwardly propagating magnetic wave which eventually dissi- 
pates in the f5'^ region. Moreover, if the oscillation frequency 
is greater than the frequency defined by equation (|36|l , part 
of the wave energy is transferred to an outwardly propagat- 
ing acoustic wave. 

Since the layer where the acoustic and the magnetic 
components are coupled is thin, and the imaginary part of 



the frequency is expected to be very small (e.g. lCunhall999l ). 
the amplitude of the oscillation changes only very slightly 
during the time that takes the wave to propagate through 
the coupling layer. Therefore, using energy conservation ex- 
pressed in terms of the dimensionless variables we can write, 
approximately. 



Fa+{ril,To) - Fa-{vl,ro) ^ Fa+{v2:ro) + F„,{7jI,to) (46) 

where to = ioi^o and io is the time at which the outwardly 
propagating wave crosses the height rj — rjl. Moreover, Fa+, 
Fa-, Fa+ and Fm are the dimensionless energy fluxes av- 
eraged over one oscillation period carried in the decoupled 
regions by the outwardly propagating acoustic component 
below the coupling layer, the inwardly propagating acoustic 
component in the same region, the outwardly propagating 
acoustic wave above the coupling layer, and the inwardly 
propagating magnetic wave in the interior, respectively. Here 
the dimensionless energy fluxes are deflned as the corre- 
sponding energy fluxes divided by the quantity R'^ujQpo. 

The average dimensionless energy flux Fa+ can be cal- 
culated using the asymptotic solution for the dimensionless 
vertic al displacement in the f3~^ region given by (e.g. iGouehl 
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Figure 6. The real parts of the dimensionless displacement perpendicular to the magnetic field (top left), along the magnetic field (top 
right), in the vertical direction (bottom left), in the horizontal direction (bottom right). Results are for a dipolar magnetic field with 
polar strength Bp = 3000 G at a latitude such that its inclination with respect to the local horizontal axis is a = 45°, and a cyclic 
frequency u = 3.08 mHz. 



(47) 



where J is the depth independent amplitude, 5 a complex 
phase and k = Rn with k the vertical wavenumber of the 
acoustic oscillation i n the interior. U s ing t his asymptotic 
form of the solution, ICunha fc Goughl (|2000l ) found for the 
outwardly propagating wave in the interior. 



\J\'^a'^ckr -^2Si -2(T,To 



8 



(48) 



where the subscripts r and i stand for the real and imaginary 
parts of the correspondent quantities. 

To determine the fraction of flux that is carried away 
by the running waves, one needs to find expressions for the 
energy fiux associated with the magnetic wave in the interior 
and the acoustic wave in the atmosphere. An exp ression for 
the former was derived bv lCunha fc Goughl (|2000D using the 
asymptotic solution for the magnetic wave given by equation 
H22p . From the latter, the authors found. 



Fm{vi,n)) 



(49) 



where D is the amplitude of the magnetic wave defined in 
equation (f22)) . 

The fiux Fa+ , on the other hand, was not considered by 



the same authors, who adopted a fully reflective boundary 
condition at the surface. To derive an expression for the 
latter, we write the average dimensionless energy flux over 
one period P carried by the outwardly propagating acoustic 
wave above the coupling layer in the form. 



2 r 



\p\v 



+ 12 



cdr 



(50) 



where is the dimensionless velocity of the oscillation, and 
the dimensionless group velocity of the wave was approxi- 
mated by c. Using equation (|4U|I . we write the solution for 
the outwardly propagating acoustic wave in the /3~ region 
in the form. 



Fi = Fu P = P II P 

Ik Ik ' 



(51) 



\v+\ = K 



dr 



r)l/2 



sin(fc|| 77 + cTrT), (52) 



where we have neglected the small imaginary part of the fre- 
quency in the amplitude, when differentiating the displace- 
ment. The corresponding dimensionless energy fiux averaged 
over one period then becomes. 



'2, TO) 



2 2 



2p 



(53) 



10 S. G. Sousa and M. S. Cunha 




0.00 



0.01 



0.02 



0.03 




-0.005 -0.004 -0.003 -0.002 -0.001 0.000 
V 




0.03 




-0.005 -0.004 -0.003 -0.002 -0.001 0.000 
V 



Figure 7. The dimensionless displacement in the interior. The 
top panel shows the vertical component while the lower panel 
shows the horizontal component. The full and dashed-dotted lines 
show, respectively, the real and imaginary parts of the displace- 
ment. Results are for the same magnetic field and the same fre- 
quency as in Fig.[6l 



Fig. [To] shows the fraction of energy flux carried by the 
running acoustic wave in the atmosphere and by the run- 
ning magnetic wave in the interior. The amplitudes \J\, \D\ 
and |A| were obtained from the numerical solutions. The 
results shown are for a dipolar magnetic field with polar 
strength Bp — 3000 G and for a cyclic oscillation frequency 
f = 3.08 mHz. Note that in this case there is a direct re- 
lation between the inclination of the magnetic field and the 
co-latitude 0, the field being in the local vertical direction 
when cos 8 = 1 and in the local horizontal direction when 
cos S = 0. As expected, the fraction of the energy fiux car- 
ried by the acoustic wave is large near the pole, where the 
magnetic field is nearly vertical. In the limit of a vertical 
magnetic field a low degree oscillation in this local analysis 
does not couple with the magnetic field. Hence, if the fre- 
quency is above the acoustic cutoff, as is the case for the ex- 
ample shown, all energy is lost. However, as the inclination of 
the magnetic field in relation to the local vertical increases, 
the fraction of energy carried away by the running acous- 
tic wave in the atmosphere diminishes, becoming negligible 
when the magnetic field becomes strongly inclined. The fact 
that no coupling takes place between the magnetic field and 
the oscillation when the magnetic field is vertical should be 
reflected also on the inwardly propagating magnetic wave 
in the interior. In fact, due to the absence of coupling no 
energy flux is expected to be carried away by this wave 
when cos 9=1. Unfortunately, numerical problems show up 
when the total fraction of energy flux carried by the running 



Figure 8. The dimensionless displacement in the atmosphere. 
The top panel shows the component perpendicular to the mag- 
netic field while the lower panel shows the component along the 
magnetic field. The full and dashed-dotted lines show, respec- 
tively, the real and imaginary parts of the displacement. Results 
are for the same magnetic field and the same frequency as in 
Fig. El 
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Figure 9. Sketch of the energy fluxes carried by the waves below 
and above the coupling region. 
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netic poles through running acoustic waves in the atmo- 
sphere and running magnetic waves in the interior. Never- 
theless, the fact that global waves with frequencies above the 
acoustic cutoff frequency are observed in some roAp stars in- 
dicates that at least part of the mode energy is conserved in 
the oscillation over each cycle. In this context it is important 
to notice that the amplitude of the magnetic component of 
the wave in the atmosphere at a given co-latitude is ap- 
proximately constant (e.g. Fig. Consequently, the total 
energy of this wave component must decrease with height, 
which means that the wave is progressively reflected. To 
check this idea we divide the magnetic wave in the atmo- 
sphere in outwardly and inwardly propagating components, 
similar to what was done for the acoustic wave in the inte- 




Figure 10. Relative energy flux carried by the running acous- 
tic wave in the atmosphere (upper panel) and by the running 
magnetic wave in the interior (lower panel), for different co- 
latitudes. The results shown are for a dipolar magnetic field with 
polar strength Bp = 3000 G and a cyclic oscillation frequency 
u = 3.08 mHz. 



waves approaches 1. Ifence we had to stop our calculations 
at cos 9 = 0.95. Still, inspection of the lower panel of Fig. [TOl 
shows that indeed the fraction of energy flux carried by the 
running magnetic wave in the interior decreases sharply as 
the fleld tends to vertical. Moreover, the energy flux carried 
away by this magnetic wave increases as the magnetic field 
becomes slightly more inclin ed, and then decreases again. 
A similar result was found bv lCunha fc Goughl l|2000'). who 
have shown that the exact position in co-latitude of the max- 
imum of energy flux carried by the running magnetic wave 
in the interior depends on the frequency of the oscillation. 

One interesting aspect to note in Fig. llUI is that the sum 
of the energy fluxes carried by the running waves is rather 
large (close to 1) for cos9 > 0.8. In practice this indicates 
that for frequencies above the acoustic cutoff, a significant 
fraction of the energy input in each cycle through the opacity 
mechanism is lost predominantly in the regions where the 
magnetic field is close to vertical. 



The expression for the outwardly propagating compo- 
nent can be derived from equation (|44|l . In the limit where 
p ^ 1 this expression becomes. 



C 

2\B\ 



(54) 



Since this wave is essentially a magnetic wave, the energy 
flux it carries while propagating outwardly averaged over 
one oscillation period is given by the time average of the 
Poynting flux. Thus, in terms of the dimensionless variables, 
we have. 



-I 



ex 



bo) 



' TO 

X (V X (3? (el^ex) x bo))], dr. (55) 



Combining equations (I54|l and (|55|l and neglecting, as be- 
fore, the small imaginary component of the frequency when 
appropriate, we find. 



, ~,2 

C a,. 



,1/2 



8\B[ 



1/2 -2<TiT-o 



(56) 

where again the amplitude |C| is to be determined by com- 
parison with the numerical solutions. From equation (|56|l 
it is immediately clear that as the outwardly running com- 
ponent of the magnetic wave travels in the atmosphere, the 
energy flux it carries decreases due to the factor p^^^. Hence, 
the total energy contained in the magnetic wave in the at- 
mosphere will also decrease with height. Fig. [ll] shows the 
fraction of energy flux carried by the outwardly component 
of the magnetic wave in the atmosphere for different heights 
in the atmosphere. Besides the rapid decline in the fraction 
of energy flux carried by this component with increasing 
height, which results from the p^^^ factor discussed above, 
it is clear that most of the energy passed onto the magnetic 
component in the atmosphere is concentrated in co-latitudes 
such that cos 9 < 0.8. 



6.2 Reflection of the magnetic wave in the 
atmosphere 

As mentioned above, for frequencies above the acoustic cut- 
off a significant fraction of energy flux is lost near the mag- 



7 DISCUSSION 

It is nowadays well accepted that pulsations in roAp stars 
are essentially acoustic in nature. Nevertheless, in the outer 
layers, where the magnetic and gas pressures are compara- 
ble, the direct effect of the magnetic field on the oscillations. 
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Figure 11. Relative energy flux carried by the outwardly prop- 
agating component of the magnetic wave in the atmosphere for 
difi'erent co-latitudes and at different heights in the atmosphere, 
namely, rj = —0.001 (full line), rj = —0.003 (long dashed line), 
7] = —0.005 (short dashed line), t] = —0.007 (dashed-dotted line). 
The results are for the same magnetic field and oscillation fre- 
quency as in Fig. 1101 



through the perturbed Lorentz force, becomes important 
and the oscillations acquire a magnetoacoustic nature. Away 
from this region, both above, where the magnetic pressure 
dominates, and below, where the gas pressure dominates, 
the magnetoacoustic wave decouples into acoustic and mag- 
netic components. Nevertheless, the presence of the coupling 
layers results in energy being transferred from the acoustic 
wave in the interior into waves that are essentially magnetic 
or acoustic in the decoupled regions, a process commonly 
known as mode conversion. 

In the present work we have used a toy model composed 
of an isothermal atmosphere matched onto a polytropic in- 
terior to investigate on the effect that mode conversion may 
have on the reflection of oscillations near the surface of roAp 
stars and also to determine the conditions under which en- 
ergy losses due to running waves are expected to be most im- 
portant. For frequencies above the acoustic cutoff, we have 
found that the maximum energy losses due to running acous- 
tic waves in the atmosphere and running magnetic waves in 
the interior takes places in regions where the magnetic field 
is not signiflcantly inclined, which for the dipolar magnetic 
field considered here, happens around the magnetic poles. 

The primary motivation of the current study was the 
understanding of the presence of oscillations with frequen- 
cies above the acoustic cutoff in some roAp stars. Hence, 
the study did not consider oscillations with frequencies be- 
low that limit. Nevertheless, from the present results and 
the results fo und by previous stud ies such as those of 
ICunha fc Gou gh (2000); Saio (2005): ICunhal ^006), which 
considered full reflection of the waves in the atmosphere, one 
can easily extrapolate the results that would be obtained 
when considering frequencies below the acoustic cutoff. In 
fact, for oscillations with frequency below the acoustic cut- 
off, the acoustic waves in the atmosphere will propagate only 
if the magnetic field is sufficiently inclined in relation to 
the local vertical. Thus, the results for oscillation with such 
frequencies are expected to resemble those found here, ex- 
cept for the fact that the fraction of energy carried away by 



the running acoustic wave in the atmosphere will decrease 
sharply to zero when the magnetic field inclination in re- 
lation to the local vertical becomes smaller than the value 
needed to ensure the propagation of the acoustic waves in 
the atmosphere. Thus, in such case, the maximal energy loss 
due to the presence of running acoustic waves should take 
place, not at the magnetic pole, but in an annulus around 
the latter, where the inner radius of that annulus is defined 
by the minimum inclination of the magnetic field needed for 
the propagation of the acoustic waves in the atmosphere. 

Despite the energy losses associated with the running 
waves discussed above, our results show that part of the 
energy is kept in the oscillation every cycle. In fact, the en- 
ergy associated with the magnetic wave in the atmosphere 
decreases significantly with height, which means that the 
wave is progressively refiected. That energy, which is kept 
in the oscillation, is concentrated in regions where the mag- 
netic field is significantly inclined in relation to the local 
vertical, i.e., around the magnetic equator. Naturally, for 
oscillations with frequencies below the acoustic cutoff, ad- 
ditional energy would be kept in the oscillation around the 
magnetic poles, where the acoustic component of the atmo- 
spheric wave would be evanescent. 

The present study considered adiabatic pulsations only, 
and it is beyond our scope to analyse the impact of the en- 
ergy losses dis cussed here on the excit ation of the oscillations 
in roAp stars. iBalmforth et al.l (|200ll ) carried out linear non- 
adiabatic calculations of pulsations in roAp stars imposing a 
transmissive boundary condition at the top when the oscil- 
lation frequencies were above the acoustic cutoff. However, 
in their calculations the authors did not consider the di- 
rect effect of the magnetic field on pulsations, which means 
they ignored both the energy associated with the magnetic 
wave in the atmosphere, which is kept in the oscillation, and 
the energy losses associated wit h the magne tic wave in the 
interior. On a more recent study. [Saicl (|2005l ) considered lin- 
ear non-adiabatic pulsations in models of roAp stars taking 
into account the direct effect of the magnetic field on the 
oscillations. However, the author imposed a fully refiective 
outer bounder condition, which means he ignored the energy 
losses associated with running acoustic waves in the atmo- 
sphere. Clearly what is needed is a study of non-adiabatic 
pulsations in models of roAp stars which takes into account 
all the aspects discussed here, including the direct effect of 
the magnetic field on the oscillations and an outer bound- 
ary condition that allows the acoustic waves to propagate 
away whenever their frequencies are above the frequency o"c 
defined in equation H36p . By providing a better insight of 
the process of mode conversion and energy exchange in a 
simple model of the outer layers of roAp stars, we hope the 
present work can serve as a motivation for the progress of 
the non-adiabatic calculations described above. 
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Figure Al. Grid of depths used in the matchings. A given value 
of the index q is assigned to a given pair of depths, such that 
when the pair of depths is in the region /3 ~ 1, g is large, and 
when the pair of depths includes a point in the f3~^ region and a 
point in the f3~ region, q is small. 



derived, respectively, from the solutions expressed by equa- 
tion (|40|) and equation ((22]), the first of which is valid in the 
P~ region and the second of which is valid in the 13'^ region. 

To compute the left hand side of the expressions given 
above, we consider the three independent solutions com- 
puted numerically, and combine them according to the rela- 
tions, 



l|3 
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(A3) 



(A4) 



From these matchings we determine the values of ci and 
C2. In practice the matchings expressed by equations ([Aip - 
(|A2|) are carried out for all depths, rather than only in the 
f3~ and (3^ regions. The depths are included in a grid with 
index q, as illustrated in Fig. lAll When /3 ~ 1 the index 
q is close to one thousand. As the index q decreases, the 
matching expressed by equation ([Aip is moved towards the 
/3~ region and, simultaneously, the matching expressed by 
equation HA2|) is moved to the /3+ region. Since the analytical 
solutions used in the matchings are valid only when q is 
small, only for relatively small values of q one may expect 
to find ci and C2 constant. These constants are then used to 
find the total solution for the displacement across the whole 
magnetic boundary layer, through the relation expressed in 
equations (|45|l . 



APPENDIX A: MATCHING PROCEDURE 

Since our calculations are linear, the amplitude of the solu- 
tion to equations p3p - (|14p cannot be determined. Thus, we 
carry out the matchings necessary to apply the outer and 
inner boundary conditions by considering the relations. 
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